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Random tensor models which display multi-critical behaviors in a remarkably simple fashion are 
presented. They come with entropy exponents 7 = (m — l)/m, similarly to multi-critical random 
branched polymers. Moreover, they are interpreted as models of hard dimers on a set of random 
lattices for the sphere in dimension three and higher. Dimers with their exclusion rules are generated 
by the different interactions between tensors, whose coupling constants are dimer activities. As an 
illustration, we describe one multi-critical point, which is interpreted as a transition between the 
dilute phase and a crystallized phase, though with negative activities. 
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INTRODUCTION 

Random matrix models [l[ provide an efficient tool to 
understand some two-dimensional statistical field theo- 
ries coupled to fluctuations of the geometry. They gen- 
erate random discretizations of surfaces. Interesting fea- 
tures we want to stress here are (i) multi-critical behav- 
iors in the thermodynamic limit, (ii) microscopic systems 
realizing them. By tuning some coupling constants, one 
achieves those multi-critical behaviors which arc differ- 
ent universality classes corresponding to conformal field 
theories (CFT) coupled to 2d quantum gravity. 

It is a remarkable fact that it is often easier to solve 
such statistical mechanical models on dynamical lattices 
than on a fixed lattice, a fact usually credited to the 
reparametrization invariance which is gained on random 
lattices (see Q for a recent application). A historical ex- 
ample is the 2d Ising model on planar graphs Q , also ex- 
actly solvable in a magnetic field [3], using a two-matrix 
model formulation which yields in the continuum the uni- 
tary Ising CFT (central charge c = 1/2) coupled to quan- 
tum gravity. 

Other multi-critical points are observed in one-matrix 
models @. The first one turns out to have the same 
entropy (or string susceptibility) exponent as the Ising 
model, which raised some confusion with the latter. But 
it was then shown that it can be derived from hard dimers 
on random lattices Q. Hard dimers are objects which 
can be attached to links of a lattice with some exclu- 
sion rule. They can be obtained from the Ising system 
at infinite temperature in an imaginary magnetic field. 
Using a matrix model formulation, the partition func- 
tion can be calculated exactly and the continuum limit is 
actually the Yang-Lee singularity, a non-unitary confor- 
mal field theory (CFT) with central charge c = —22/5 
coupled to two-dimensional quantum gravity (see also 
@; HI)- That universality class is found generically in 
models of hard objects [9(. We want to emphasize that 
the derivation from a microscopic system (hard dimers) 



helped identify the continuum limit. The mechanism is 
that dimers generate an effective action for a one-matrix 
model with a new interaction potential. Hence, that pro- 
vided generic one-matrix potentials with a physical pic- 
ture, corresponding to non-unitary matter CFTs coupled 
to quantum gravity [H, [j| EH 11 1 . 



To generate random lattices in dimensions d > 2, one 
can use random rank-c? tensor models 13; 14 L A class 
of such models (known as colored models 15-3]) is now 
known to have a continuum limit in the large N scaling 
limit (N being the size of each tensor entry) 1 191 . While 
different universality classes have been found [20|, which 
are multi-critical points in the matrix model terminol- 
ogy, they (i) come from a toy model which reduces to a 
rectangular-matrix model, (ii) have not been physically 
motivated so far. 

First guesses could be to run the Ising model or hard 
dimers, on those random lattices. However unlike in two 
dimensions, using the large N limit of a colored tensor 
model it has been found that there is no phase transition 
at finite temperature 21[ in agreement with previous nu- 
merical analysis [22j |. 

In this article, we introduce a family of colored ten- 
sor models which exhibit multi-critical behaviors in a re- 
markably simple fashion. A way to characterize them 
is the entropy exponent, which quantifies the prolifera- 
tion of microscopic configurations in the thermodynamic 
limit. We find the same entropy exponents as those of 
[20I ] which also corresponds to the entropy exponents of 
multi-critical branched polymers [HI ■ 

Moreover, we provide a physical picture of such mod- 
els. The presence of different interactions between ten- 
sors is interpreted in terms of dimcr insertions on the 
dominant graphs, with some generalized exclusion rules. 
The coupling constants are the dimer activities. In the 
dilute phase, one recovers the entropy exponent of pure 
random lattices, strengthening its universality. Like in 
matrix models, multi-critical behaviors appear when ac- 
tivities are appropriately tuned. However, the thermo- 
dynamic limit is defined in a slightly different way from 
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the standard one where it is strictly the number of lattice 
sites which is sent to infinity. We show that some multi- 
critical points correspond to phase transitions from the 
dilute phase to a crystallized phase, though with nega- 
tive activities. Finally, we discuss other interpretations 
of our model. 




RANDOM TENSORS AND HARD DIMERS 

Random tensors and random melonic triangulations 
of the sphere 

We start with the simplest colored tensor model 
which generates random triangulations of d-dimensional 
pseudo-manifolds. We take d + 1 pairs of rank-d tensors 
( T ni...n d ^n 1 ... nd )c=o,...,d, each carrying a color index c, 
and having components n< = 1, . . . , N. The tensor ac- 
tion has a quadratic part, 



Slink — - ^ 

9 ^ 



E 

c=o {m=i,...,N, 
i=l,...,d\ 



1 n 1 ...n d ± n\...n d i 



(1) 



and an interaction part with a positive vertex and nega- 
tive vertex, 
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The coupling constant g will be used to reach the contin- 
uum limit. The free energy 
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can be written as a sum over closed, connected, (d + 1) 
colored graphs, 



m = E E 



j^-ui(G 2n ) gdn 



(4) 



n£N colored graphs 

G271 with In nodes 



where m(G) is a p ositive integer known as the degree of 
the graph [la4lq . Those colored graphs are built from 
(d+l)-valent nodes of positive and negative orientations, 
respectively generated by S+ and 5_ . It means that the 
half-links around each node carry the colors to d, and 
are drawn on the plane by going clockwise or counter- 
clockwise. The propagator coming from S^nk can only 
join a positive node to negative one, by connecting two 
half-links of the same color. 

Those graphs encode triangulations of rf-dimcnsional 
pscudo manifolds in the following way. First, we iden- 
tify fc-cclls, k = 0, . . . , d, on the graph as the connected, 



Figure 1. On the left, the elementary melon with external 
color 0. On the right a typical melonic graph, obtained by 
inserting melons inside melons. 



closed sub-graphs with exactly k colors |l5(. Then, it is 
found that there is a colored triangulation which is dual 
to each graph, obtained by replacing each (d — k)-ce\\ 
with a fc-simplex. For example, a node is dual to a d- 
simplex and a link between two nodes corresponds to a 
(d — l)-simplcx shared by two neighboring d-simplices. 
The boundary maps of the triangulations arc defined by 
the boundary maps of the colored graph HBj . 

In the large N limit defined in 16l4l8|. only triangu- 
lations of the d-sphere contribute. Moreover, they must 
have ui(G) = 0, from which they inherit specific combi- 
natorics described in [l9|, where they have been coined 
melonic graphs. The basic building block is the elemen- 
tary melon, formed by two nodes of opposite orientations 
which are joined by d links, as shown in the figure[T] The 
two external lines have the same color. The graphs are 
formed by gluing series of elementary melons next to each 
other, and inserting such sequences inside one another. 
They are melons inside melons so to say. To get a closed 
graph (hence a triangulation of the <i-sphere) we join the 
two external lines. A typical graph is depicted in the 
figure m and we refer to [19| for details. 

As a gluing of simplices, the elementary melon is built 
by gluing two d-simplices according to the rules, along 
d of their boundary simplices. It has two external {d — 
l)-simplices of the same color which corresponds to the 
two external lines of the dual graph. If one thinks of 
the simplices as being flat and equilateral, then melonic 
triangulations carry strong positive curvatures. 

The continuum limit is the regime where F is domi- 
nated by graphs (triangulations) with a very large num- 
ber of nodes (d-simplices). It is reached by sending the 
parameter q to a critical value g c , the radius of conver- 
gence of F [19(. There F becomes singular and behaves 
like F ~ (g c - 5 ) 2 "T with 7 = 1/2. 



Hard dimers 

A dimer is an object which can be attached to a link 
of a lattice. We consider hard dimers, with a general- 
ized exclusion principle. It is usually considered that two 
dimers can never touch each other. Here we choose that 
a node can share dimers with only one neighboring node. 
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Figure 2. A graph dual to a triangulation at d = 3 (since 
nodes are four-valent) with a set of two dimers on the left, 
a single dimer in the center and a set of three dimers on the 
right. 
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Figure 3. The dimer parts of the action can be drawn with 
internal (dashed) lines which represent dimers and internal 
faces (2 lines and 1 face for the set of two dimers on the left, 
3 lines and 3 faces for the set of three dimers on the right). 
The effective interaction only sees the external legs carrying 
some T and T. 



It is reasonable if we think of a dimer as a bond between 
two monomers. If two nodes with monomers are joined by 
several links, they can form a bond on each link, provided 
they do not share a dimer with another partner. Inter- 
preting further, two dimers joining the same two nodes 
define a surface whose boundary is formed by the dimers 
and along which the dimers may couple to each other. If 
a third dimer is added, one gets three surfaces which glue 
together to form a closed surface. It is bounding a 3-cell 
which can be interpreted as a three-dimensional coupling 
region. Such coupling regions can have dimension d at 
most. An example is given in the figured! An activity 
Zfe is assigned to each set of k dimers joining two nodes, 
and the partition function is 

E Gn ( Zk ,g)= £ g-^Uz^. (5) 

DeV(G n ) k 

Here D is a dimer configuration, i.e. an assignment of 
dimers to the lattice satisfying our exclusion principle, 
and T>(G n ) is the set of dimer configurations. \D^\ is the 
number of sets of k dimers joining two nodes in D and 
we have set \D\ = J^k ^|-Dfc| the total number of dimers 
in D. At this stage g is just a book-keeping parameter 
which can be absorbed into a redefinition of the activities. 

Making the lattice dynamical means that we sum over 
a set {G n } of graphs. If colored graphs are considered, 
that can be implemented using the above tensor model 
supplemented with effective interactions which create 
hard dimers. Links without dimers are created by Sunk, 
positive and negative nodes without dimers by S± . A set 
of d+l — k dimers with colors {0, . . . , d}\{a, . . . , Ck) be- 
tween two nodes is obtained by connecting a vertex + and 
a vertex — with links of colors {0, . . . ,d}\ {ci, . . . , c>}. 
The external lines have colors {ci, . . . , c&} and carry in- 
dices determined by the form of 5 + , S- . The scaling with 
N is also determined by S+, S- and by the ( d + 1 - k )( d - k ) 
internal faces formed by the dimers. That leads to effec- 
tive interactions, i.e. from the view of the external lines 
which do not carry dimers, of the form 

fe 

n{ci,...,c k } _ Zd+l-k T~\ rpa^iCi (( ,\ 

'-'dimers — (k-l)(2d-k) 11 X X 1 1 ■ V s ) 

9 N 2 . =1 

We have denoted the index contractions with a trace, 
as it is otherwise a bit cumbersome. We have depicted 



•^dimers and ^dimers for d = 3 on tne n g ure El The terms 
with k = 1 are just corrections to Sunk- 
The free energy S(zfc,<7) is defined by 

e JT^(» l , t ) = t Y[dT*dre- s , 

-° P) 

S = Sunk + S+ + S- + ^ X! '^di'm'crs^'' • 

fe=l {ci,...,Cfc} 

For simplicity, we set z\ = 0. Sdirners does not affect 
melonicity of the large N limit but instead really intro- 
duces dimers on melons. A typical Feynman graph is a 
closed connected melon Gin with a dimer configuration 
D and the free energy has the expansion 

E(z k ,g)=J29 dn E S G 2 „(^,.9). (8) 

nGN {G 2 „} 

It is interpreted as a grand-canonical partition function 
for dimers, where one sums over all possible melonic tri- 
angulations and all dimer configurations. There is an as- 
sociated canonical partition function where the quantity 
p = dn — \D\ is kept fixed, 

Sp(*) = E E E IK* 1 , w 

new {G 2re } £>er>(G 2 „), fc 

\D\=dn-p 

such that 

^,9)^/S p (4 (10) 

The corresponding thermodynamic limit is the large p 
regime. Note that it is not the standard regime of lat- 
tices with an infinite number of sites (studied in Q for 
instance), as p = dn — \D\ also depends on the number 
of dimers. Equivalent ly, In g is not a chemical potential 
for the number of lattices sites n, but for p instead. The 



1 It is possible to solve the model similarly with Z\ ^ 0. However, 
this becomes more complicated and does not bring different qual- 
itative results. 
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large n regime will be studied in a companion paper, to 
appear. 

To make sense, E p (zk) should behave like 



E p (z k ) ~ Ap^ 3 g c (z k ) p , 



(11) 



for some constant A and some function g c {zk)- The expo- 
nential decay is expected since {g c {zk)\ is then the radius 
of convergence of E(zk,g). For a fixed set of activities 
{zfc}, when g is tuned to g c (zk), the grand-canonical par- 
tition function E,(z,g) starts to be dominated by lattices 
with very large values of p (again, this is slightly different 
from being dominated by lattices with a very large num- 
ber of sites). The power-law decay p 7-3 then determines 
the singular behavior of the system close to that critical 
surface, 



3(zk,fl) ~ (.9 - 9c{zk)) 2 1 - 



(12) 



7 is known as the entropy exponent, as it is related to 
the way the microscopic configurations involved in S p 
proliferate at large p (also known in two dimensions as 
the string susceptibility exponent, for which one-matrix 
models give 7 = — 1/m). 



THERMODYNAMIC LIMIT AND 
MULTI-CRITICAL BEHAVIORS 



A key object of the analysis is the full, connected 2- 
point function which we write 



(TZ ... nd %... n , d )=U(z k ,g) Ip„ 



(13) 



Then, 
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In the melonic sector, all the above expectation values 
are monomials in U. Part of the proof can be found in 
24|, and more details will be appear soon elsewhere. In 
11], it was shown that (S+) ~ (S-) ~ U d+1 . More- 
over, (^aimere' *^) is obtained by reconnecting the exter- 
nal T Ci with the external T Ci using the full 2-point func- 
tion, which gives a U k behavior. Hence if F has a singular 
part like (|12[) , the most singular terms in the above equa- 
tion are those going like U, with U ~ (g c (zk) — g) 1 ^ 1 ■ 

U can be evaluated as the geometric series of the 1- 
particle irreducible 2-point function S, U = g/(l — gT,), 
say with external color 0. A self-consistency equation 
for U is obtained if one can express £([/) independently. 
It is easy in the melonic sector. One contribution comes 
from the dimcr-frcc nodes, S±. We glue a vertex + with a 




Figure 4. In the melonic sector the 1PI 2-point function is 
built from the full 2-point function U and dimer insertions. 
The first term of the right hand side is that for pure random 
lattices [ljj, due to S±, and the others are due to effective 
sets of dimers between the two external lines. 



vertex — by inserting d full 2-point functions, one on each 
color c = l,...,d. The terms Sdimera a l so contribute. 
One isolates the two lines of color as the external lines, 
and connect the other external lines which have the same 
color inserting U . That gives 



S = -\L 
9 



k=0 



(15) 



also symbolically depicted in the figure 2] for d — 3. The 
closed equation on U and g is 



9 = (1 - z d )U 



d-2 

E 

k=l 



(16) 



By definition, U is a series in g and the thermodynamic 
limit corresponds to the regime where g goes to the ra- 
dius of convergence of U. The criterion for that loss of 
analyticity is 



93 
dU 



0. 



(17) 



which is polynomial in U of degree d (in particular, it has 
explicit solutions for d = 3,4). When the activities are 
small enough, the physical root is that passing throug 
U c {z k = 0) = l/(d+ lf/ d , and satisfying 



d ~ 2 /rl\ 
l-z d = (d+ l)U d + ]T z d _ k (k + 1) I I) U 

k=l ^ ' 



(18) 



That gives a critical surface z d (z2, ■ ■ ■ , z d —i, U). Start- 
ing from the physical solution U(zk = 0), it can be in- 
verted to give U(zk), hence g(z k ), as long as dzd/dll 7^ 0. 
Along this surface 7 = 1/2 and it shows the universality 
of the critical behavior of pure random melonic lattices 
observed in (l9j . 

The condition dzd/dU = is 

- 2dz d _ 1 = (d+ ^dU' 1 - 1 + z d-k(k + U k -\ 

k=2 ^ ' 

(19) 



2 This is the model studied in figjl though with different conven- 
tions. 
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which restricted to the critical surface ([18)) gives one re- 
lation between the activities. Since dzd/dU = is equiv- 
alent to d 2 g/dll 2 — 0, one has (U — Uq) ~ [g — g c (zk)) ^ 
and hence 7 = 2/3 on this sub-surface. 

It actually holds as long as dz d _i/dU ^ 0. When the 
derivative is zero, we get a second relation between the 
activities which corresponds to d 3 g/dU 3 = 0, so that 
7 = 3/4. This reasoning can be repeated so that at each 
step the activities Zd, ■ ■ ■ , z m can be expressed in terms 

Of Zm-l, ■ ■■ ,Z 2 . 

As dg/dU is a polynomial, that set of relations means 
that some of its roots meet. When the activities are var- 
ied, the physical solution may intersect with other roots 
and jump from the universality class of pure random lat- 
tices to a new one. If exactly m = 2, ...,d — 1 roots 
meet, 



d k g 
dU k 



= 0, k = 1, 



,m, 



but 



i) 



rn+l 



dV 



r ^ °- (2°) 



Those relations determine m — 1 activities as functions of 
the others. On such points U ~ (<? — <7 C ) m+1 ^! so that 
multi-critical points with 



1 



7 



for to = 2, . . . , d, 



(21) 



are obtained. Like in two dimensions [6j], that requires 
some activities to be negative, hence "non-physical" . 

The behavior of the system may be complicated in the 
neighborhood of such points. There are two notable situ- 
ations we want to mention, (i) The cut singularity, when 
it is not possible to go beyond some values of the pa- 
rameters. A well-known example in two dimensions is 
the Yang-Lee singularity which is indeed the tri-critical 
point of the matrix model for hard dimers Q, and also 
appears generically in models of hard objects @. Tech- 
nically, it is due to the collapse of the physical branch 
with another root of dg/dU such that both solutions be- 
come complex beyond the singularity, (ii) The phase 
transition, when the physical branch meet another root 
of dg/dU and that they exchange their role after the crit- 
ical point. A famous example is the 2d Ising transition 
on random lattices at finite temperature observed in a 
two-matrix model Q and also in models of hard objects 
on bicolorable random lattices Q. 



PHASE TRANSITIONS 



natural to go to higher dimension^. In the remaining we 
consider d = 6 with z\ = Z3 = Z5 = 0, so that 



*L = _ 7f/ 6 

du 



75z 2 C/ 4 -45z 4 C/ 2 + l-z 6 . 



(22) 



For small activities, it is a phase of random lattices with 
dilute dimers. Therefore it can be described perturba- 
tively, and U grows linearly with the activities. For large 
activities, we expect a phase of fully packed dimers. The 
main contributions then come from Sdimers, which means 
that the term —7U e , coming from nodes without dimers, 
can be neglected. Assuming the dimer terms all have 
the same order, one finds typically that U/y/z~2 goes to 
zero. It is actually expected that the free energy does 
not change too much with zu in the crystallized phase. 

When all activities are zero, U = (1/7) 1 / 6 . First, we 
lower zq from to —6, so that U = 1. The equation (fT6|) 
for g(U) has two real extrema, the physical continuum 
limit being the positive one, and correspondingly dg/dU 
has to real roots, ±1, and four complex roots e 2lk ^^, 
k = 1,2,4,5. To observe multi-criticality, it is necessary 
to create two real positive roots, and to collapse one of 
them with the physical root. Therefore we tune some 
activities to create two real positive roots, setting z 2 = 
— 1| , Z4 = y| and z 6 = —111 for instance. 

Then we find a path to collapse and exchange the role 
of the physical root with another one. For example, 



91 5 

za{z 2 ) = "60 + 28 Z 2(25z 2 - 14), 
375 

z 6 (z 2 ) = -174 - — z 2 (15z 2 + 14). 
The real positive roots on this path arc 



„ , , 1 3 / 50 
U x {z 2 ) = 2, 



. , 1 3/50 



(23) 
(24) 

(25) 
(26) 

(27) 



For -4^ < z 2 < — |i, one has U Q < U\ < U 2 and the 
physical solution is Uq. However, Uo and U\ meet when 
z 2 reaches the critical value z 2 = and 7 = 2/3 at 

this point. For z 2 < z 2 , U\ < U < U 2 . 

To see there is a phase transition, notice that Uq ^ 
y/—Z2 at large negative activities, so that it is not a 
physically reasonable solution. By contrast, Uf <C z 2 Uf, 



Note that the coefficients of the polynomial dg /dU can 
be chosen arbitrarily except for the vanishing of the term 
U^ 1 which corresponds to the constraint that the sum 
of the roots is zero. That leaves enough freedom to ob- 
serve the typical behaviors mentioned above. While it is 
possible to do it with d = 3 or 4, it is easier and more 



To observe a phase transition with d = 4, one has to increase 24 
above 1, which means that the quadratic term in the action has a 
vanishing coupling at some point and then comes with the wrong 
sign. While it seems there is no trouble with the self-consistency 
equation, we do not know to which extent going through 24 = 1 
is problematic. 
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(a) U(z2) with a discontinuity 
in 42- at z% 



(b) U(z2) with a discontinuity 

d 3 U , * 

in -jTT at z 2 



Figure 5. Plots of U{z2) for different paths to the phase tran- 
sition. Before the transition, Z2 > z 2 , the physical solution is 
Uo (normal style), and for Z2 < z 2 , it is U\ (dashed). 



which means that the constant solution is meaningful 
when nodes with dimers dominate. We conclude that U\ 
is the solution for z 2 < z 2 and that there is a crystalliza- 
tion transition between a dilute phase and a phase with 
fully packed dimers. 

In the continuum limit the free energy is a function of 
Z2, F(z2, 24(22), zq(z2), (7(22)). The order of the transi- 
tion is determined by the discontinuity in some deriva- 
tive of F with respect to 22- The same way dF/dg can 
be written in term of U, dF j dz-i is also a polynomial in 
U. According to the path we have chosen, U is obvi- 



ously continuous at the transition (see figure 5(a) I, but 



its derivative is not, and therefore the transition is second 
order. 

However, it is possible to choose other paths in the 
space of activities to make the transition as smooth as one 
wants. Let us give an example. Take a parametrization 
such that Uq = 2 + x for zi > z%, the transition being at 
x = 0. Instead of setting U\ as a constant like above, one 
can make it behave like U\ ~ 2+x — x 3 around x = 0. For 
definiteness, choose U2 = 3 + x and the negative roots as 
the opposite of the positive ones. One can expand such a 
polynomial to identify Z2{x), £4 (a;), zq(x). Z2 is of order 6 
in x but we can use some truncation to extract x(z2) in 
the neighborhood of z 2 = —119/75. Since the change of 
variables between x and z 2 is well-defined, one concludes 
that only d 3 U/dz 2 is discontinuous at the transition, so 
that the transition is of the fourth order (see figure [5(b)] ). 



APPLICATIONS 
Monomers on random branched polymers 

The system of dimers on random melonic triangula- 
tions can be mapped to a problem of monomers on ran- 
dom branched polymers. That is suggested by the self- 
consistency equation (fl"6|) which is similar to the alge- 
braic e quat ions on trees which are also found in matrix 
models 



polymers, due to a partial order on the set of melons in a 
melonic graph [l9| . Those branched polymers are rooted 
(d+ l)-ary trees, as any 1PI connected 2-point subgraph 
of the graph dual to a triangulation is represented by a 
vertex of degree (d + 2). A vertex of the tree either is a 
leaf (it has no child) or it has d + 1 descendants. 

Dimers with their exclusion rules are simply mapped to 
monomers attached to branched polymer vertices which 
are not leaves. We have d—1 types of monomers, coming 
with activities^ Z2, ■ . ■ , Z4. While the use of monomers on 
branched polymers trivializes the dimer exclusion rules, 
there is one subtlety which actually comes from the col- 
oring. Remember that the dual graph to a triangulation 
has colors on its links. They become colors on the links 
of the corresponding tree. A vertex comes from its par- 
ent with a link of color, say, i, and its descendants are 
joined with links of colors to d, including the color i. 
The rule is that such a vertex can carry one monomer 
with activity z 2 , Z3, . . . or Zk provided it has at least k 
children which are leaves with colors different of i. 

Intuitively, one can think of such monomers as being 
able to grasp some leaves attached to a vertex, but being 
sensitive to the color of their parent, they cannot grasp 
a leaf with the same color. 



Two-dimensional interpretation 

Finally, we would like to draw attention towards a dif- 
ferent interpretation of our model. The triangulations we 
built are obviously <i-dimcnsional. However, their dual 
graphs in the melonic sector are planar graphs. More 
precisely, they form a subset of the planar graphs en- 
countered in matrix models. This suggests a interpreta- 
tion on the 2-sphere, but we are not able to describe it 
yet. Remarkably, going to this subset of planar graphs 
changes the universality classes observed in 2d gravity 
coupled to matter. The 2d entropy exponents are of the 
form = — 1/m (with m = 2 for pure random lat- 
tices, i.e. pure gravity) while we got families of the form 
7 = l+72d (771 = 2 being again pure random lattices fl9j|). 
Technically, one can count planar graphs with two exter- 
nal le gs n ot necessarily in the same face using decorated 
trees 



llj . resulting in the same asymptotic behavior as 



for melons. But the free energy is not directly related to 
it by a simple derivative. 

It is known that matrix models at multi-critical points 
correspond to conformal field theories coupled to 2d 
quantum gravity, and it would be interesting to find 
whether the same statistical models on melons also take 
part to the conformal field theory landscape. A first hint 
would be to rewrite the Schwinger-Dyson equations de- 
rived in [25j as Virasoro constraints. 
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First, there is a one-to-one correspondence 
between melonic triangulations and a family of branched 



We again assume for simplicity z\ = 0. 
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